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Abstract 

The ARGO-YBJ experiment has been in stable data taking from November 2007 till February 
2013 at the YangBaJing Cosmic Ray Observatory (4300 m a.s.L). The detector consists of a sin¬ 
gle layer of Resistive Plate Chambers (RPCs) (6700 m^) operated in streamer mode. The signal 
pick-up is obtained by means of strips facing one side of the gas volume. The digital readout of 
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the signals, while allows a high space-time resolution in the shower front reconstruction, limits 
the measurable energy to a few hundred TeV. In order to fully investigate the 1-10 PeV region, an 
analog readout has been implemented by instrumenting each RPC with two large size electrodes 
facing the other side of the gas volume. Since December 2009 the RPC charge readout has been 
in operation on the entire central carpet (~ 5800 m^). In this configuration the detector is able to 
measure the particle density at the core position where it ranges from tens to many thousands of 
particles per m^. Thus ARGO-YBJ provides a highly detailed image of the charge component at 
the core of air showers. In this paper we describe the analog readout of RPCs in ARGO-YBJ and 
discuss both the performance of the system and the physical impact on the EAS measurements. 

Keywords: Air shower detection, RPC detector. Calorimetry 


1. Introduction 

The ARGO-YBJ experiment, promoted and funded by the Italian Institute for Nuclear Physics 
(INFN) and by the Chinese Academy of Sciences (CAS) in the framework of the ItalyChina sci¬ 
entific cooperation, has operated at the YangBaJing Cosmic Ray Laboratory (Tibet, PR China, 
4300 m a.s.l., 606 g/cm^) from November 2007 till February 2013. As typical of an air shower 
array apparatus, it has benefited from a large field of view (~ 2 sr) and a high duty cycle (> 85%), 
that allowed a continuous monitoring of the sky in the declination band from -10° to 70°. The two 
key features, namely the dense sampling active area and the operation at high altitude, allowed 
the study of the cosmic radiation at an energy threshold of a few hundred GeV. 

ARGO-YBJ has been proposed as an experiment capable of investigating a wide range of 
fundamental issues in Cosmic Ray and Astroparticle Physics at a relatively low energy threshold; 
high energy y-ray astronomy, at an energy threshold of a few hundred GeV; search for emission 
of y-ray bursts in the full GeV-TeV energy range; Cosmic Ray (CR) physics (energy spectrum, 
chemical composition, pjp ratio measurement, shower space-time structure, multicore events, p- 
air and pp cross section measurement) starting from TeV energies; Sun and heliosphere physics 
above 1 GeV. 

The detector could measure in different energy ranges, according to the way of operation, 
which essentially means the minimum measurable spatial density of particles. In scaler mode, in 
order to lower the energy threshold down to 1 GeV m , the total counts on sub-units of 44 m^ are 
measured every 0.5 s, with limited information on both space distribution and arrival direction 
of the detected particles. In shower mode, position and arrival time of the charged particles 
are measured, then the arrival direction of the primary particle is determined by a fast timing 
technique [|2]- 

In shower mode ARGO-YBJ can determine the number of particles in the shower either by 
counting the number of fired pick-up strips (digital readout) or measuring the total charge induced 
by the particles passing through the detector (analog readout). The maximum particle density 
which can be measured is different for the two cases, simply depending on geometry in the first 
case (23 strips/m^) or on electronics in the second case (up to many 10"^ /m^), corresponding to 
a maximum detectable energy of a few hundred TeV in the first case and to many PeV in the 
second case. 
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From 2007 to the end of operation ARGO-YBJ has been operated in scaler mode and in 
shower mode. At December 2009 the analog readout was made operational so increasing the 
maximum achievable energy up to many PeV. Due to the dimension of the detector, the gamma 
astronomy studies are limited to energies < 20 TeV 0, well addressed by the digital readout 
(see i) and reference therein). Therefore, in what follows, we will consider the major results 
obtained so far by ARGO-YBJ on the CR physics and the new opportunities offered by the analog 
mode of detector operation. 

• Using the strip multiplicity spectrum measured by ARGO-YBJ the light-component (pro¬ 
ton and helium) spectrum of the primary CR in the energy region 5-200 TeV has been 
evaluated by using a Bayesian approach 0, so bridging for the first time direct and indi¬ 
rect measurements below 100 TeV, region not accessible by other EAS experiments. The 
ARGO-YBJ data agree remarkably well (within about 15%) with the values obtained by 
adding up the proton and helium fluxes measured by CREAM Ena, both concerning the 
total intensities and the spectrum slope, so confirming the CREAM finding that the proton 
and helium spectra, from 2.5 to 250 TeV, are both flatter compared with the lower energy 
measurements. Later on the energy spectrum of proton and helium has been measured be¬ 
low the so-called ’’knee”, from 100 to 700 TeV, by using a hybrid experiment|8] made of 
a wide fleld-of-view Cherenkov telescope 19] and the ARGO-YBJ carpet. The light com¬ 
ponent has been well separated from other CR components by using a multi-parameter 
technique where the analog information has been used for the first time. A highly uniform 
energy resolution of about 25% was achieved throughout the whole energy range (100- 
700 TeV). The measured spectrum agrees in both spectral index and absolute flux with the 
spectrum obtained by ARGO-YBJ alone in the lower energy range from 5 TeV to 200 TeV. 
Since proton and helium nuclei are the bulk of the CR at energies below the knee (3x lO'^ 
eV), the study of their spectrum in this energy region is of primary importance. 

• Basing on the CR flux attenuation for different atmospheric depths, i.e. zenith angles, and 
exploiting the high detector accuracy in reconstructing the shower properties at ground, 
ARGO-YBJ has measured the production cross section between 1-100 TeV CR protons 
and Air nuclei ITOl . The analysis results have been also used to estimate the total proton- 
proton cross section at center-of-mass energies between 70 and 500 GeV. The total proton- 
proton cross section has then been inferred from the measured proton-Air production cross 
section by using the Glauber theory. The result is consistent with the general trend of 
experimental data, favoring an asymptotic In^(s) rise of the cross section. 

• ARGO-YBJ has reported IHTI the observation of anisotropic structures on a medium an¬ 
gular scale as wide as ~ 10° - 45°, in the energy range 10'^ - 10^^ eV. The intensity spans 
from 10 to 10“^, depending on the selected energy interval and sky region. For the first 
time, the observation of new MSA structures throughout the right ascension region 195° - 
290° has been observed with a statistical significance above 5 s.d.. At higher energies both 
the EAS-TOP IfT^ and IceCube lfT3l experiments observed significant anisotropy around 
400 TeV. At this energy, the signal looks quite different from the modulation observed up 
to 50 TeV, both in amplitude and phase; it suggests that the global anisotropy may be the 
superposition of different contributions from phenomena at different distances from the 
Earth. 

• In the ARGO-YBJ experiment, for the first time, the geomagnetic held (GeoME) lateral 
stretching is observed M at small distances from the shower core. This effect has been 
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suggest by Cocconi ns who supposed that the lateral displacement induced by the Earths 
magnetic field is not negligible with respect to the Coulomb scattering when the shower 
is young. The non-uniformity in the azimuthal distribution has been deeply studied by 
ARGO-YBJ concluding that it is well described by two harmonics: the first one of about 
1.5% and the second one of about 0.5%. The first harmonic is due to the GeoMF; the 
second one is the sum of the magnetic and detector effects. The phase of the first harmonic 
can be used as a marker of the absolute pointing of an EAS array and, if any, can be used to 
apply a simple correction to the absolute pointing. The ARGO-YBJ experimental results 
represent the first evidence of the charge density reduction near the EAS axis due to the 
GeoME. 

With the analog readout in operation on 5800 m^ ARGO-YBJ can extend its CR physics 
results at higher energies (> 100 TeV), so reaching the knee region of the CR spectrum, and cope 
with new physics items. The energy spectrum of all CR primaries, the all-particle spectrum, 
will be measured up to the knee region and the evolution of the different components will be 
investigated too with special attention to the light component. In this respect we notice the 
methodological difference with respect to the sampling apparatus, in fact ARGO-YBJ will use 
only information provided by the very near region of the core, that essentially means truncated 
size, local age and particle density at the core position. 

Improvements are expected on the p-Air cross section measurement, that is the extension to 
center-of-mass energies in the TeV region. The possibility of performing this measurement actu¬ 
ally strongly relies on the knowledge of the chemical composition around the knee, specifically 
of the light component, which is still matter of debate and investigations. So the light compo¬ 
nent and the p-Air cross section measurements look as intimately correlated and will represent 
a key item in the ARGO-YBJ investigations. Considering the geometry of the detector, the core 
region is measured with unprecedented detail and the very forward region of the interaction can 
be studied. It will be interesting to see how well the different models will describe the data at 
higher energies. In this context it is worthwhile to note that the proton-air inelastic cross section 
measured by the ARGO-YBJ experiment in the energy range 1-100 TeV has been found in good 
agreement with the values set in the CORSIKA/QGSJET code, and that, according to the results 
shown in||5l, the QGSJET and SYBILL models provide the same description of the longitudinal 
development of the shower below a few hundred TeV. 

Due to the size of the detector, then to the statistics of collected data, it will be almost impos¬ 
sible to detect the CR anisotropy signal at high energies especially considering the level of 10 
reported ifT^ . even using not perfectly internal events. 

The study of the GeoMF effect on the shower development are expected to be performed at 
higher energies, in fact an azimuthal modulation was already observed not only at the Yakutsk 
array for EAS with energies above 50 PeV UB and at the ALBORZ Observatory for energies 
above 100 TeV IfTTl . but also in radio experiments IfTSl . 

Shower events with multiple core structure, or multicore events, have been observed and 
investigated by many experiments with different techniques, namely emulsion chambers EH 
and hadronic calorimeter EOll . Those events are typically explained within the framework of jet 
production, which is essentially provided by the leading particle interaction with the Air target 
nuclei, the separation between the cores being related to the pr of the leading particle. The 
ARGO-YBJ carpet with its huge surface (~ 5800m^) allows the measurement of events with 
cores separated up to 100 m, so providing a good opportunity to measure very high pj events 
ED. Events at high pj are a perfect tool to investigate hadronic interaction models or even new 
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physics. In this respect also the study of double front showers, at high energy, would be a check 
for some exotic physics Il22ll . Measuring the shower core region will allow the study of local 
clustering phenomena of particles in the shower that could have a strong correlation with the 
mass of the primary hadron; in fact this phenomena should be enhanced around the core region 
ll23l . In this paper the detector layout is described in §2 while the main features of the analog 
signal from the RPC are discussed in §3. The §4 copes with the question of the intrinsic linearity 
of the RPC detector and results are reported on a specific test-beam. Then in §5 the trigger of 
the experiment and the related electronics are introduced. The architecture of the charge readout 
system and the local trigger that enables the charge readout on the single Cluster are discussed 
in §6, while in §7 the logic operations of the charge readout system are presented focusing on 
synchronization, data collection and data packing. The main operations of the DAQ control 
system concerning the analog readout (managing, calibration and monitoring) are described in 
§8. Then, in §9, the calibration chain of the analog information is presented and the stability of 
the detector is discussed. Checks of consistency and the final performance of the analog readout 
of the detector are shown in §10. A summary of the experimental results is given in §11 while 
the conclusions are drawn in §12. 

2. The ARGO-YBJ detector 


O 
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Central Carpet 130 Clusters 


Guard Ring 23 Clusters 



(1.39 X 1.23 m2) 


Figure 1: The ARGO-YBJ detector setup. The Cluster (12 RPCs) is the basic module. The Local Station (LS) is the 
DAQ unit providing local trigger and readout. 


The ARGO-YBJ detector ll24ll . hosted in a building at the YangBaJing Cosmic Ray Observa¬ 
tory, consists of a central carpet ~ 78x74 m^, made of a single layer of Resistive Plate Chambers 
(RPCs) with ~93% active coverage, enclosed by a guard ring partially instrumented (~20%) up 
to ~110x100 m^. The apparatus has a modular structure, the basic data acquisition unit being a 
Cluster (7.6x 5.7 m^), made of 12 RPCs (1.25 x2.80m^ each). The full detector has 153 Clusters 
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(130 in the central carpet, 23 in the guard ring) with a total active surface of about 6700 m^, as 
shown in Fi^ Each RPC is read out by means of 80 pick-up strips (61.8x6.75 cm^, the spatial 
pixels) facing one side of the gas volume. The fast-OR signal of 8 contiguous strips defines the 
logical pad (61.8 x 55.6 cm the time pixel) which is used for timing and triggering purposes. 
Any manifold coincidence of fired pads of the central carpet (Npad) above a given multiplicity 
threshold (Ntrig), namely Npad ^ Ntrig in a time window of 420 ns, implements the inclusive trig¬ 
ger that starts the event data acquisition. The apparatus, in its full configuration of 153 Clusters, 
has been in smooth and stable data taking since November 2007 till February 2013 with a duty 
cycle > 86%. The trigger threshold is Npad ^ 20 (Eq > 300(900) GeV for primary photons (pro¬ 
tons)) and the corresponding trigger rate is ~ 3.5 kHz. The high granularity of the detector and 
its time resolution provide a detailed three-dimensional reconstruction of the shower front. The 
digital pick-up of the RPC, which has a density of 23 strips/m^, can be used to study the primary 
spectrum up to energies of a few hundred TeV; above these energies its response saturates ll25l . 
In order to extend the measurable energy range and fully investigate PeV energies, where particle 
densities at the core position are larger than 10^/m^, each RPC has also been equipped with two 
large size electrodes of dimensions 1.23 x 1.39 m^. These pick-up electrodes, called Big Pads 
(BP), face the other side of the RPC gas volume (see Fig|^ and provide a signal whose amplitude 
is expected to be proportional to the number of charged particles impinging on the detector. 

The BP and the 40 corresponding strips see the same signals provided by particles impinging 
on the RPC and detected by it. The RPCs are operated in streamer mode, with a gas mixture of 
Argon (15%), Isobutane (10%) and Tetraflouroethane (75%). The operating voltage is 7.2 kV. 
This setting provides a typical efficiency > 95% with an intrinsic time resolution of about 1 ns 
and a minimum ionizing particle (m.i.p.) signal on the BP of ~2 mV. 


3. Chamber Layout and Signal Characteristics 

The chamber layout and the results concerning the digital read-out have already been dis¬ 
cussed in Il24l . According to the layout shown in Figj^the BP is a large pick up electrode 
grounded through a 1 kfl resistor in order to avoid slow charging, typical of a floating electrode, 
followed by a breakdown discharge with possible electronics damage. 

It is made by a copper foil (17 fim thick) glued on a PET foil (190 fim). 

A typical signal of the BP is shown in Figj^ it has a rise time in the range of a few tens of 
ns and a discharge time of about 5 yus; the amplitude of the signal, on a 50 fl load, ranges from 
millivolts to tens of Volt. In Fig|^the amplitude distribution of the signals generated by a m.i.p., 
as measured at the experiment site, is reported along with the gaussian fits to the first and second 
peak of the distribution. The ’’single particle” trigger was provided by a small telescope f2S\ 
made of three RPCs of 50x50 cm^, two above and one below the chamber under measurement, 
with about 1 cm of lead put on top of the lower RPC, and the triple coincidence of their signals 
used to select a m.i.p.. The amplitude distribution of the signal provided by the BP on a 50 Q 
load, measured by a ChargeMeter board specifically designed to this aim, shows a first peak 
at ~ 2.5 mV which represents about 90% of the sample and corresponds to one streamer, with 
cr ~ 0.4 mV, and a second peak at about 4.4 mV with cr ~ 0.6 mV) , which corresponds to 
two streamers. The first populated bin shows a spike at 0.7 mV which corresponds to electronic 
noise, in fact its content is consistent with the overall RPC inefficiency, ~ 3%. As reported in 
Il24l . a foam layer of a few millimiters has been added between the RPC and the BP in order to 
avoid the coupling with the strip plane mounted on the other side of the RPC. In Fi^is reported 
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Figure 2: Layout of the Argo-YBJ RPC with High Voltage, Low Voltage and grounding connections (see for details). 


the mean amplitude of the m.i.p. signal versus the applied HV, with and without the foam layer. 
It can be seen that in case of no foam there is a much faster increase of the measured amplitude, 
which is higher by a factor 2-^3 in the efficiency plateau; the foam instead stabilizes the signal 
and makes the plateau extension more evident. Similar behavior has been observed for some 
parameters like strip multiplicity and pulse height of the strips with respect to the thickness of 
the foam layer Il24l . 

4. Intrinsic linearity of the RPC 

In order to check the intrinsic linearity of the RPC response up to particle densities of a few 
10"^/m^, which is the density at the core of showers induced by primary cosmic rays of PeV 
energies, a test-beam on small RPCs has been performed in October 2009 at the Beam Test 
Facility (BTFIEtI. part of the DAFNE-factory complex (INFN National Laboratory of Frascati, 
Italy). The facility includes a high current electron and positron LINAC, a 510 MeV e“ and e^ 
accumulator and two 510 MeV storage rings. The BTF is a beam transfer line optimized for 
the production of electron or positron bunches, in a wide range of multiplicities down to single¬ 
electron mode, in the energy range between 50 and 800 MeV. The typical pulse duration is 10 ns, 
very close to the time thickness of the shower front near the core, and the maximum repetition 
rate is 50 Hz. 

The experimental setup of the test at BTF is shown in Fig|^ Four RPCs of dimensions 
63 X 57 cm^, having the same layout as the standard ones (see Il24l and Fig|^, in particular 
provided with 8 strips for digital readout on one side of the gas volume and one pad for ana¬ 
log readout on the other side, were put just at the exit of the beam line. With all transfer line 
quadrupoles off during the test, the beam spot was completely de-focused and limited, at the exit, 
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Figure 3: Typical signal from the Big Pad on a 50 H load. The signal is reported in a) with a horizontal scale of 2.5 fis 
and in b) with a horizontal scale of 100 ns. In both cases the vertical scale is 10 mV. In Fig[^ the vertical scale is 
200 mV (horizontal scale 2.5 fis) with the oscilloscope set at infinite persistence and triggered for 20 minutes by particles 
crossing the chamber. 


by the vacuum pipe section which has dimensions 5x3 cm^. The bunch rate was 1 Hz and the 
slits defining the beam characteristics were operated so to have a number of particles ranging 
from a few up to many tens per bunch. The energy was in the range 450-500 MeV. The RPCs 
were operated with the same mixture used in the ARGO-YBJ experiment and a check of their 
performance was done before the test beam by stacking the 4 RPCs in a vertical telescope, then 
measuring the efficiency and the operation voltage with cosmic muons crossing the telescope. 
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Figure 4: Amplitude distribution of the Big Pad signals induced by m.i.p., at 7.2 kV, as measured at the experimental site 
of YangBaJing . The pedestal (the BP electronic noise) as well as the single and double streamer peaks are clearly shown. 
The bin size is ~0.1 mV.The two curves represent the best gaussian fits to the first and second peak of the distribution. 


Above the knee voltage, ~ 9.2 kV, in the plateau region, the efficiency was > 95%. Both values 
were in agreement with the expectations. At 9.2 kV, the voltage of operation during the test, the 
measured amplitude of the pad signal for a single particle crossing the detector was 3.7 + 0.1 
mV. A small lead glass calorimeter, located just beyond the RPCs (see Fig|^ and aligned with 
its major axis along the beam direction, was used to measure the particles exiting the RPCs. The 
lead glass block is from the former OPAL experiment Il28l and has the shape of a truncated prism 
of Schott SF57 lead glass; the block length is 37 cm, corresponding to 24 radiation lengths, with 
a base of 11 x 11 cm^. The light signal is read by a Hamamatsu R2238 photomultiplier. The data 
acquisition was performed by a custom system managed by a Motorola 6100 CPU and housed 
in a VME crate. The calorimeter signal was used for cross comparison with the RPC signals, 
particularly with the signal of the hrst RPC along the beam line. It was acquired, after suitable 
adaptation, by means of the same electronics used for the RPC analog signals. The trigger to the 
DAQ was provided by the signal of bunch-crossing, available through the BTF system. The beam 
geometry (refer to m for details), even with the beam fully de-focused, guaranteed an almost 
complete containment of the beam by the calorimeter, or very marginal loss of information (the 
Moliere radius of the Schott SF57 lead glass is 2.61 cm). The scatter plot of Fig|^ shows the 
signal of the hrst RPC along the beam line (mean value + r.m.s.) in bins of ADC counts in the 
calorimeter. Assuming the linearity of the RPC behavior and the amplitude of the single particle 
mentioned before, we reconstructed between 7 to 30 particles impinging on the RPC surface. 
This number is fully consistent with the estimate provided by the beam monitoring system of 
the BTF. On the right scale of Fig|7]the corresponding particle density of the beam is reported. 
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Figure 5: Mean Amplitude versus applied HV, with/without foam, as measured at the experimental site of YangBaJing. 


To check the consistency with linear response, the experimental data have been fitted with the 
red straight line shown in Fig|7]and the residual values, normalized to the fit values, are reported 
in the histogram of Figj^ The gaussian fit to the residual distribution (Figj^ shows a good 
agreement, as confirmed by the value of the Id-o.f.. From the fitted values of the gaussian 
parameters one can say that local deviations are contained within a few per cent (r.m.s) , while 
the integral deviation (mean) is below 1%. The offset of the RPC response in Fig|^is due to 
the strong attenuation of the calorimeter signal and to its adaptation to match the specifications 
of the readout electronics. In conclusion, up to 30 particles on 15 cm^ there is no evidence of 
deviation from linearity behavior of the RPC, which means linearity response up to density of 
about 2 X 10'*/m^. Of course this value is conservative because the particle density of the beam 
spot is not properly uniform. 

5. Local Station and Trigger System 

The trigger of the experiment is generated by the digital signals sent by the Front-End boards 
mounted on the RPCs. These digital signals are processed by a specific crate named Local Station 
(LS) Il29ll - the Cluster DAQ Unit -, as depicted in Fig. that provides the pad multiplicity to 
the trigger system. The LS crate contains and manages 12 receiver cards, one I/O card for the 
communication with the DAQ and one active backplane. Each receiver card collects the signals 
coming from one RPC chamber and provides the fast-OR signals which start the TDC counting. 
When a trigger occurs, a common stop signal goes from the backplane to the receiver cards, 
which store the patterns of the active strips, and to the TDCs which perform the arrival time 
measurement. Each LS outputs two busses, namely a 6-bit Low Multiplicity (LM) weighted bus 
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Figure 6: Experimental setup of the test beam at the BTF. 


(providing signal when >1, >2, >3, >4, >5, >6 pads are fired within 150 ns) and a 4-bit High 
Multiplicity (HM) weighted bus (providing signal when >7, >16, >32, >64 pads are bred within 
60 ns). 

The trigger system |[30l, which has the LM and HM busses in input, implements two different 
selection algorithms based on a simple, yet robust, majority logic which takes into account the 
topology and the time distribution of the bred pads. The LM trigger implements a selection of 
small-size showers by requiring at least Ntrig bred pads in the central carpet. The HM bigger 
has been designed to select showers with a much higher particle density. The whole bigger 
electronics is hosted in 2 VME crates as shown in Fig. 

The LM bigger is implemented with a four-level hierarchical architecture, where each level 
correlates only pads belonging to adjacent areas. According to simulations and to the measured 
pad rate (~400 Hz/pad), the number of spurious hits in the 420 ns trigger window has been 
estimated to be less than 3. 

The data collected in each LS, that is the pattern of the bred ships and the arrival time of the 
particles, are packed and transmitted at each trigger occurrence to the central DAQ at a rate of 
160 Mbit/s (16-bit word in 100 ns) by means of the I/O card. 

The present trigger set-up enables just the LM selection with a threshold of 20 pads. Since the 
amount of data for each event strongly depends on the shower size and the cosmic ray spectrum 
follows a power law, the data frame of the event ranges from about hundred bytes to Mbytes, 
with an average event size of about 2 kbytes. 

6. The Charge Readout System 

The BP signals of two adjacent Clusters are processed by elecbonic modules hosted in a cus¬ 
tom crate, called MINICRATE, that has two independent sections, each one containing 3 readout 
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Figure 7: Result of the RPC linearity test performed at the BTF (see text for details). The fit with a straight line, in red, 
has been performed. 


reg. value (Rfs) 

7 

6 

5 

4 

3 

2 

1 

0 

f.s. voltage (V) 

0.29 

0.58 

1.16 

2.13 

4.26 

8.52 

16.2 

32.4 


Table 1: 3-bit register values and corresponding effective full scale (f.s.) voltages; each f.s. voltage has a 
spread of 2-r2.5% due to the different board productions. The scale of operation is indicated by Gn with n 
equal to the value of R^. 


cards (CHargeMeter cards) and a Control Module (see FigfT0)i. The CHargeMeter (CHM) pro¬ 
cesses 8 analog signals and digitizes them, while the Control Module builds the data frame of 
3 CHM boards and transfers it to the LS, which finally provides the data to the central DAQ. A 
12.5 m coaxial cable is used to feed the signal to the CHM input adapted to 50 Q. The CHM 
board is made of 8 identical sections. The first block of each section consists of a voltage limiter 
able to cut out spikes greater than 40 V to protect the next electronic stages. It is followed by 
a linear amplifier with a dynamic gain set by a 3-bit programmable register, which provides 8 
input values for full scale (f.s.) setting, as show in Tabjl] Operating the system at different scales 
allows an overlap between digital and analog readout which is an efficient way to calibrate the 
analog system. The scale of operation is indicated by Gn with n equal to the value of Rfs as 
shown in Tab[T] 

The sensitivity of each scale is about 0.2% of the full scale value, below this value the electronics 
sensitivity is reduced. 

The signal at the output of the amplifier ranges from 0 to 2.5 V which is the reference input of 
the Analog to Digital Converter (ADC). 

In the shaper module the long tail of the BP signal is cut out by doubling and inverting the signal, 
then summing the original to the inverted one with a At delay. This time delay has been set to 
about 500 ns mainly to guarantee that the analog signal reaches its maximum and, secondly, to 
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Figure 8: Residuals normalized to the fit values (see Fig|2). The gaussian fit to the distribution shows a good agreement 
(see text for details). 



Figure 9: Schematic drawing of the Front-End electronics and DAQ system. The FE block is based on a 8-channel 
discriminator chip for strip signal processing located inside the metallic casing of the RPC. The MINICRATE block is 
in charge of the BP signal processing. The MINICRATE electronics and the LS crate are located in the middle of the 
Cluster unit; each LS in turn is connected to the Central Station with a star-like custom network for triggering and data 
transfer purposes. 


be sensitive to delayed particles in the shower front. A Peak and Hold (P&H) circuit represents 
the core of the next stage. It continuously samples the output of the previous stage and keeps the 
highest reached value for 2 yus; the voltage drop per time unit of the P&H is ~ 4 mV/jUs, which 
does not cause any trouble with respect to the trigger jitter time as already shown in ||3T|| . If a 
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Figure 10: The logic block diagram of a MINICRATE section. 


conversion signal (named local trigger) arrives within 2 jt/s, the ADC starts to digitize the P&.H 
amplitude. The conversion signal is generated by the HM bus of each LS and can be selected 
with different thresholds (namely >16, > 32 , > 64 and > 73 fired pads). 

If no local trigger arrives within 2 ^s, the P&.H is reset and then it starts again to sample the 
signal from the previous stage. 

The ADC data collection is managed in each section by the Control Module via a custom bus 
protocol operated on the backplane lines. The Control Module receives the local trigger from the 
LS and distributes it to the 3 CHM units starting the ADC digitization. Then, each CHM board 
replies by asserting a local Busy signal in order to prevent the generation of further local triggers 
while the ADC is converting. After 14 clock cycles at 10 MHz the data are converted, latched 
and ready to be collected by the Control Module. 

In this project the AD7472 component ll^ by Analog Devices has been used. It is a 12 bit 
high-speed ADC offering an high throughput with a low power consumption (1.5 MSPS with 4 
mW power consumption). The differential nonlinearity (DNL) measured by the producer is 1.9 
LSB. Some tests about performance, homogeneity of the board production, homogeneity of the 
electronic channels, linearity and calibration features have already been done and results have 
been reported in IISTI . ||33], Il34| . ll^ and ll^ . 


7. The Data Flow 

In this section we describe the analog data flow, from the FE electronics up to the event 
builder, paying attention to the synchronization with the digital frame and the overall system. 
After the digitization of the BP signals, the data are sent to the FIFO of a special receiver board 
(Reel 1 in Fig. hosted in the LS. The transfer to the central DAQ starts when the experiment 
trigger confirms the local trigger. A multiplexer logic provides a simple and flexible data readout. 
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Figure 11: Average values of Total Dead Time (TDT) and Analog Dead Time (ADT) versus the local trigger rate. 


Each CHM board readout is performed in 4 clock cycles (400 ns, thus the whole information is 
transferred in 1.2 fis to the FIFO of Reel 1 in the FS. If the local trigger is not validated by the 
experiment trigger, the data stored in Reel 1 are discarded and the ADC system will be ready to 
process a new local trigger. In both the cases the Busy signals generated by CHM boards are 
reset. 

When the experiment trigger arrives at the FS, the event building process starts and the event 
number, the addresses of the fired strips, the pad time information and the BP charge information 
are packed in a data frame which is transmitted to the DAQ system, then stored into the FIFO 
memory of the Data Buffer board llJTl of the Fevel-1 (FI) readout system ( IIMl , 1391). A simpli¬ 
fied block diagram of the data flow is shown in Fig. 


The ARGO-YBJ DAQ modular structure allows a high speed and efficient data collection. 
It is built on a two-layer read-out architecture implementing an event-driven data collection 
by using two custom bus protocols, based on VME-bus HOl . A Motorola VME Processor 
MVME6100, is responsible for the complete read-out of a Fevel 2 chain. Three MVME6100 
CPUs are used to manage and control the trigger and the charge readout system. 

In the last part of this section we show the result of a test on the data collection performance with 
different analog configurations. Fig. 11 shows the behaviours of the Analog Dead Time (ADT) 
and the Total Dead Time (TDT) as a function of the local trigger rate. The fraction of the TDT 
due to the ADT is less than 10% for a local trigger rate of IkHz, then it increases up to 60% at 
4 kHz. As the local trigger rate for the ARGO-YBJ detector does not exceed 100 Hz, it is clear 
that the analog readout does not introduce any significant dead time; we estimate the ADT to be 
less than 2% of the TDT at the experimental conditions. 
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8. The Run Control 


The Run Control has been developed in order to provide the interface between the shift oper¬ 
ator and the DAQ components for managing and monitoring the run. The Run Control is based 
on the following modules: the run manager which allows the operator to startup a new run; the 
log of the inter-processes communicator that provides messages coming from the DAQ nodes; 
the run monitoring that shows some run information like run number, trigger issued, time left, 
dead time etc. The Run Control system has been successfully used to control the DAQ system. 
The commands pertaining to the charge readout system, issued by the Run Control both for hard¬ 
ware setup and monitoring/debugging tasks, are distributed to the Analog VME CPU as shown 
in Fig. 1^ All the MINICRATEs are connected in a daisy-chain network to the controller boards 
hosted in the VME crate allowing a two-way serial communication line at the speed of 1 Mbit/s. 
The processes pertaining to the CHM boards are managed by the Control Module which sets the 
ADC full scale of operation and the local multiplicity threshold for triggering. 

Since the CHM boards, located near each Cluster, are distributed over a very wide area they could 
suffer from variations of the environmental parameters, therefore their calibration and continuous 
control are fundamental. The calibration runs are started by the Run Control: a full scale is set, 
then signals with known amplitude are generated and sent to the ADC inputs, finally the digitized 
values are collected and flagged properly. A Digital-to-Analog Converter (DAC902 HTl l from 
Burr-Brown Products by Texas Instruments, hosted in the Control Module, generates pulses of 
variable amplitudes. When a calibration pulse is generated, a fake HM signal is also generated in 
the LS at same time. This signal is processed by the trigger system which provides the Common 
Stop signal that starts the readout on the CHM board. The DAC902 is a high-speed DAC offer¬ 
ing a 12-bit resolution. To improve the signal/noise ratio, the differential configuration output 
in conjunction with the operational amplifier has been used. The characterization of the DAC 
device shows a linear response in the whole range of 12 bits. Its typical differential nonlinearity 
(DNL) is 0.5 LSB. 


9. Calibration and Stability of the Detector 

In this section we discuss both calibration and stability of the RPC analog detector. The de¬ 
tector performance is shown, the steps of the calibration procedure are described and the related 
results discussed. 

9.1. The analog data taking 

The analog system has been completely assembled at YangBaJing in its final layout on the 
130 central Clusters and then put in data taking since December 2009. 

The system has been operated at f.s. 0.29 V (Rfs=7, maximum gain) till June 2010. At 
this scale the sampled particle density overlaps with the particle density measured by the digital 
readout system. These data have been used both to study the detector behavior and for calibration 
purposes. From July to middle August 2010, the system was operated at the intermediate scale 
corresponding to f.s. 2.13 V (Rfs=4); since middle August 2010 the f.s. 16.2 V (Rfs=l) was set. 
The number of BPs in the central carpet is 3120, being 3% the number of dead or badly working 
channels. A typical air shower event with core hitting the detector is shown in Fig[^ In the 
vertical scale is reported the number of particles measured in each BigPad, whose area is about 
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Figure 12: Shower as recorded by ARGO-YBJ through the analog readout system. According to Montecarlo simulations, 
the shower would correspond to a primary proton of about 2 PeV. In the vertical scale is reported the particle density 
measured in the BigPad, whose area is about 1.7 m^; nx and ny are the coordinate of the BigPad. The area with low 
particle counts coiTesponds to the central carpet(~78x74 m^.) 


1.7 m^. According to Montecarlo simulations for proton primary, the shower has an energy of ~ 
2 PeV. 

The lateral distribution function of this shower is reported in Fig|T^ Each blue point is the 
measured density of particles on the specific BP located at distance r from the core, which has 
been calculated as the center of gravity of the BP signal distribution. The black points and the 
their error bars correspond to mean value + r.m.s. of the particle density measured in bins of 
distance from the shower core. The figure shows the performance of the analog RPC readout: 
the study of the core region of high energy showers can reach an unprecedented richness of 
details. Owing to the analog readout, the measurement of the proton-air cross section IfTOll and 
the composition studies Q will be extended to PeV energies. 

9.2. The gas distribution and the BP signals 

Up to the installation of the analog system in the ARGO-YBJ experiment, RPCs have been 
mainly used for triggering and tracking purposes. On the other hand, gas detectors operated in 
quasi-proportional mode ll42ll or in streamer mode ||43]| have been successfully used in calorime¬ 
ters. This is the first time RPCs are used in analog readout mode over a very large area. The gas 
distribution is accomplished by feeding parallel lines (’gas channels’) each one made of a pair 
of RPCs serially connected, as shown in Fig[^ The 4 BPs along the gas line are labeled with 
BP0,BP1,BP2 and BPS. A total of 780 gas channels supply the RPCs carpet while additional 
lines provide the gas distribution to the guard-ring RPCs. Each gas line is equipped in input with 
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Figure 13: Lateral distribution of the particle density for the event shown in Fig[T^ The red point are the measured 
density of particles on the BPs, while the black points are the mean values (± r.m.s.) of the particle density measured in 
bins of distance, 0.5 m wide, from the shower core. 


a capillary tube (glass, 32 mm long and 0.6 mm diameter) whose impedance is so high to prevail 
on every other one and to define by itself the gas flow through the gas channel. The gas volume of 
the 153 Clusters is 19 m^, and the gas flow is set to 3.1 volume/day; the effective gas flow, evalu¬ 
ated through the gas consumption, or by measuring the rate of weight decrease of the gas bottles, 
corresponds to 3.0 volume/day. An effect related to the gas flow has been observed, namely 
the amplitude decrease of the BP signal along the gas line. This effect, which is particularly 
strong between BPO and BPl, is shown in Fig[T^where the amplitude distributions are reported 
separately for the BPO and BP3. The unit of measurement of the amplitude, mV/strip, assumes 
one fired strip to correspond to one detected particle, this point will be discusses and clarified 
later on (Sec. VIII.III.2). The maximum attenuation (~17% ) is observed for BP3 with respect to 
BPO. Moreover, while the mean amplitudes of BPl, BP2 and BP3 are tightly correlated to the 
environmental parameters, specifically pressure and temperature, the amplitude of BPO spans in 
a band wider than the environmental parameters allow M- This effect was unexpected and is 
still matter of investigation. 

9.3. The calibration procedure 

In order to convert the ADC counts to particle number two steps are needed, namely first 
converting the ADC count to amplitude, then the amplitude to the number of particles. The first 
step is referred to as electronic calibration, the second one as gain calibration. 
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Figure 14: Gas distribution: the elemental ”gas channel” is shown. It consists of two RPCs serially fluxed. Inside the 
RPCs the gas flow is represented by the thick red arrow; the gas inlet is split in two pipes which feed the flrst chamber 
(lower left), then two pipes connect the exits of the first to the inputs of the second chamber, and finally the two output 
pipes are joint to form the gas outlet. Along the gas line the 4 BPs are labeled as BPO, BPl, BP2 and BP3, BPO being the 
first pad and BP3 the last one receiving the gas. 


9.3.1. Electronic calibration 

Calibration runs llMl l45l have been performed at the different f.s. the electronics has been 
operated; the relation between ADC count and input amplitude, provided by the DAC, has been 
fitted with a suited functional form as shown in Fig|^ The coefficient of the linear term (PI), 
which is the dominant one, comes out to have a spread 1.5-2.5% depending on the board pro¬ 
duction, which in fact was achieved in three different bunches. In Figj^the percent error of the 
calibration is reported versus the ADC count for the 3 scales of operation. At the f.s. 16.2 V 
and 2.13 V, the error is by far below 1% for ADC count greater than 100 while it is less than 
2-3% below 100; at the f.s. .29 V the error is below 1% for ADC count greater than 400 and 
becomes 5-6% below 100, in between it is 2-3%. To check the long term stability of the electron¬ 
ics, calibration runs have been issued several times after weeks and months; the results showed a 
stability within about 1 % . 


9.3.2. Gain calibration 

The gain calibration has been performed by using the data taken with f s. .29 V, which is the 
most sensitive one and allows a direct comparison between analog and digital measurements of 
the same quantity, that is the number of particles hitting the RPC. Vertical showers {6 < 15°) 
have been selected. The adopted calibration procedure relies on the number of fired strips as a 
good estimate of the number of particles crossing the detector. The good linearity between the 
number of fired strips and the number of detected particles has been verified by simple simulation 
of the detector response. To include the signal induction on adjacent strips (as expected in case 
of particle crossing nearby the strip edge) a strip multiplicity 1.10+0.02 per particle, as directly 
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Figure 15: Gain (mV/strip), or amplitude of the BP signal con'esponding to one fired strip. The gain distribution for the 
l.st (BPO, red) and the 4.th (Bp3, blue) BP along the gas line is reported. 


measured at the experiment site, has been assumed in the simulation. The Fig[T^ shows the 
dependence of the number of bred strips on the number of detected particles in the range of 
interest to the gain calibration, the shown band corresponds to a variation of the strip multiplicity 
between 1.08 and 1.12. The typical plot of the BP signal amplitude (Amp) versus the number 
of fired strips is presented in Fig|T^ a good linearity is shown up to 15-18 fired strips, while a 
clear deviation can be seen even below 7-8 bred strips. Accordingly a linear bt was performed 
in the range 8 to 15 strips that maximize the overlap region with linear response, so debning the 
BP gain (G(mV/strip)) as the slope of the btting line; 

C(„.V/«„» = (1) 

fired-strip 

in case of dead ships in the BP area the bt range has been suitably re-scaled. This btting proce¬ 
dure has been applied to calibrate each BP. The number of collected events is enough to reach a 
statistical uncertainty in the gain determination lower than 2% That can be achieved with about 
1 hour of data taking. 

Then, taking into account the relation between bred strips and detected particles, shown in 
Fig|^ the gain G can be redebned and expressed in mV/particle; 

Cl.,,.VIpanicle) = X ( 2 ) 

ired-Strip detected-particle 

In the following, since the relation shown in FigfT^applies to all BPs, G is typically reported 
in mV/strip. 
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Figure 16: Electronics calibration of a channel. The relation between ADC counts (output) and signal amplitude (input) 
of the specific channel has been obtained by fitting the experimental points (stars), provided by the calibration system 
(see text), with a suited functional form. 


The gain distributions of BPO and BP3 have been reported in Figjl^ Dividing the data sample 
in 1 hour sub-samples and repeating the procedure all over the sub-samples, the gain variation 
with time has been obtained. The dependence on atmospheric pressure (P) and temperature (T) 
has been investigated. These environmental parameters are measured with an high accuracy, 
respectively of +0.5mbar and +0.25K. In order to display the gain dependence on P and T, the 
mean gain has been calculated for each group of BPs having the same position along the gas 
line. As an example, the behavior of the mean gain of BPS, along 8 contiguous days, is shown in 
Fig|^ As expected the gain is daily modulated and correlated to both atmospheric pressure and 
temperature. 

In Figj^the best fit to the data it is shown for the relation 

G(0 = G(Oo X ^ X ^ (3) 

where Pq and To are arbitrarily chosen and G(Oo is the gain of specific BP at the reference Pq 
and Tq, which has been used to correct the data for changes of the environmental parameters T 
and P. The same degree of correlation has been observed for BPl and BP2, while BPO appears to 
suffer of wider variations. 

By correcting the gain of each BP for the dependence on P and T, the spread of the gain 
distribution clearly reduces and we are left with a spread of about 2%. The gain variation after 
correction has been used to evaluate the stability of the detector over long time. In Fig|^ the 
BP3 gain is reported versus time. The horizontal scale has a break corresponding to about 2 
years. The points before the break refer to about three months of data, while the last two point 
refer to about 3 days of data taken after 2 years. The variation corresponds to a gain increase of 
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Figure 17: The percent difference of the functional form with respect to the calibration points is reported versus the ADC 
count for 3 scales of operation. 


+2.5%/year; same values have been found for BPl and BP2 while for BPO we have +9.5%/year, 
therefore a mean derivative of +4.4%/year due to the higher variation of BPO. 

As seen, at the most sensitive scale (Rfs=l, f.s. 0.29 V) this procedure guarantees a direct 
way to calibrate the analog system providing to each BP channel its own gain, on a run by run 
base, owing to the overlap between analog and digital information. 


9.3.3. All Scales Calibration 

At higher f.s., the procedure is quite similar; first the electronic calibration is performed, 
then the gain measured by the 0.29 V scale is used to extract the number of particles. The gain is 
taken at Pq and Tq which correspond to the mean values of the environmental parameters in the 
specific ran and the correction for P and T variations is applied. This procedure relies on detector 
stability, which is certainly true on a single run base. The long term variations are taken into 
account directly by the gain modulation as obtained by the calibration procedure, while for very 
long times also the gain derivative with time, shown in Fig|^ is considered. All information 
needed to extract the number of particles (n) hitting the specific BP, that is electronics calibration 
parameters and gain of each channel, along with the correction for pressure and temperature, is 
managed through a Data Base. Given the amplitude A of the measured signal (in mV) and the 
gain G (mV/part), the number of particles hitting the specific BP is simply obtained as n - AjG 
with an error that, assuming A and G do combine linearly, can be written as: 


O'n _ crA CTq 

n A G 


(4) 


According to the described procedure and to the equation Q, many sources contribute to the 
error on G, namely: 
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Figure 18: Relation between the number of fired strips and the number of detected particles, with (blue squares) and 
without (red squares) induction on adjacent strips as it results from Montecarlo simulations. Blue squares refer to strip 
multiplicity 1.1 (that is 10% probability of signal induction on adjacent strip), the band (dashed lines) corresponds to a 
variation of the strip multiplicity between 1.08 and 1.12 (± 1 s.d.). 


1. the error in the fitting procedure for the mV/strip determination 

2. the error coming from the strip versus particle correlation function (e^p); 

3. the error on the electronic calibration (eecai)', 

4. the error on the P/T correction {epTcor)- 

The specihc errors for each source, both statistical and systematic, are conveniently summarized 
in Tab|^ Therefore for the gain G of the BP we have (JcIG - (2.3%)stat + (3.8%)^^^. 

Summing in quadrature statistical and systematic errors we get 4.4% as total error on the 
gain. As far the amplitude, given the r.m.s. of the single m.i.p. signal (Sec. II) and the linearity 
of the RPC response (Sec. Ill), we have 


(Ta 20% 

A -sfri 

provided the measurement error on A is small, like in our case. Finally we have 


n 


20% 

yfn 


+ 4.4% 


(5) 

(6) 


The dependence of the gain on the zenith angle of the detected air showers has also been 
investigated. Data reported in Fig|^ show that there is no dependence on the incidence angle. 
Actually, given the width of the gas gap, 2 mm, and the typical size of the streamer (~3 mm), 
one can expect almost no dependence at least up to 50-60 degrees. 
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Figure 19: BP amplitude versus number of fired strips. To obtain the gain a linear fit has been performed in the range 
8-15 strips (heavy line). The slope of the fitting line defines the gain G(mV/strip). 


10. Performance Evaluation. 

In order to estimate the goodness of the calibration procedure the data concerning each BP 
have been analyzed at the different steps of the calibration procedure, in terms of ADC count, 
amplitude and number of particle. The distribution of each of these quantities, on a single BP, 
follows a power law as a consequence of the shape of the primary cosmic ray spectrum. There¬ 
fore, by htting the distribution, the power index of the mentioned variables has been measured 
for each BP. If the calibration procedure operates correctly, then the spread of the index distri¬ 
bution reduces going from the raw ADC count to the reconstructed number of particles. This 
check has been done for both f.s. 0.29 V and f.s 16.3 V, the results are summarized in Tabj^ 
Looking at values of the table a few considerations can be done, namely; a) the spread of the 



Table 2: Error sources (%) to the determination of G 


24 














Figure 20: The BP3 gain (mV/strip) monitored throughout 8 days; temperature and pressure, suitably rescaled (a.u.) to 
fit the gain scale, are also reported. 


slope distributions in case of the ADC count are almost independent of f.s. (ADC count row); b) 
the electronic calibration (Ampl row) seems to perform better at higher f.s. but this is due just to 
the calibration procedure, in fact more points are used at higher f.s.; c) in case of strips (particles) 
the slopes and their r.m.s. values are quite independent of the scale. Finally, according to the last 
row, 4.3% is a good estimate of the homogeneity among BPs after calibration. 

In order to check any systematic effects related to the use of different scales, which means the 
correct matching between different energy ranges, events with core in a fiducial area of the cen¬ 
tral carpet (the 6x9 central Clusters, corresponding to 2380 m^) and zenith angle 6 < 15° have 
been selected and the distribution of the number of particles on the BP at the shower core position 
(PMax), or the BP with the highest signal, has been studied. PMax is correlated to both energy 
and mass of the primary particle. The differential rate of Pmax is reported in Figj^for the two 


<> ±rms 

(nnsi <>) 

f.s. .29 V 

f.s. 16.2 V 

ADCcou 

Ampl 

Part 

2.39 ±0.19 (8.2%) 
2.55 ±0.18 (7.1%) 
2.51 ±0.11 (4.4%) 

1.90 ±0.18 (9.4%) 
2.56 ±0.13 (5.1%) 
2.51 ±0.11 (4.2%) 


Table 3: Index of the power law distribution of ADC count (ADCcou), Amplitude (Ampl) and number of 
particles (Part) in a BP: The mean value and the r.m.s. of the index distribution (see text for details), for all 
BPs, in case of f.s. 0.29 V and f.s. 16.2 V, are reported. 
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Figure 21: The BP3 gain versus the T/P parameter. A linear relation fits the data quite well, it is used to correct the data 
for changes of the environmental parameters T and P. 


scales of operation G4 and Gl. It can be seen a perfect matching between contiguous ranges and 
absence of any systematic effect. We notice that Pmax spans over two and half decades, while 
the event frequency runs over five decades. 

An absolute comparison of the experimental distribution in Figj^has been done with Monte- 
carlo expectations. For this aim, showers have been generated by the Corsika ll46l code with 
QGSJETII as hadronic interaction model. Energies have been generated according to the spec¬ 
tra in HU for proton, and up to Ee, primaries. The induced showers have been sampled at the 
YangBaJing elevation. The shower cores were distributed in an area wider than the fiducial area 
of the central carpet (64 x 64 m^) and almost vertical showers have been considered {9 < 30°). 
Detector and trigger have been fully simulated. The same cuts used for the real event selection 
were applied to the Montecarlo sample. The Montecarlo expectation for PMax is reported in 
Eigj^for an absolute comparison with data, showing a very good agreement. The same events 
have been used to estimate the energy threshold of the analog trigger with respect to different 
primaries. The trigger efficiency versus energy, for different primaries, namely protons. He and 
Ee is reported in Eig|^ it shows that above 100 TeV the analog system is fully efficient for all 
primary particles. The relation of PMax to the energy of the primary has been obtained by the 
above mentioned Montecarlo simulations at fixed energies, the results are reported in Eigj^for 
protons. He and Ee. Considering the RPC intrinsic linearity, which has been demonstrated at 
least up to particle density of 2x10"^ m“^, and the electronics features, this implies we can fully 
explore the CR spectrum up to ~ 10 PeV. The related statistics of detectable events, per year of 
operation and vertical events, is shown in Tab|^for different primaries . 
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Figure 22: BP3 gain behavior versus time. The horizontal scale has a break coiTesponding to about 2 years. The points 
before the break refer to about three months of data, while the last two point refer to about 3 days of data taken after 2 
years. The gain variation corresponds to a derivative of the detector response of +2.5%/year. 


Primary 

E,a (TevT’-'--^^ 

P 

He 

Fe 

100 

3.08e+05 

2.89e+05 

1.35e-t05 

500 

1.81e+04 

2.02e+04 

1.04e-t04 

1000 

4.82e+03 

6.07e+03 

3.45e-t03 

2000 

1.07e+03 

1.68e+03 

1.14e-t03 

4000 

1.70e-H02 

3.91e+02 

3.75e+02 

10000 

8.5 

3.37e-t01 

8.2e-t-01 


Table 4: Number of vertical {9 < 15 deg) events per year of operation in a fiducial area of ~ 2400 for 
E > E,i, and different primaries. 


11. Summary 

In this paper the analog detector of the ARGO-YBJ experiment and its performance have 
been presented. The features of the RPC signals are discussed within the framework of the DAQ 
of the experiment. Results on the intrinsic linearity of the RPCs, as measured at the BTF of 
the INFN National Laboratory in Frascati (Italy), are reported. These results confirm a linear 
behavior up to particle density of 2 x 10"^ nT^, at least. The analog detector is able to operate at 
different energy ranges, owing to the flexibility of the readout electronics. At the most sensitive 
scale, the digital readout and the analog readout perfectly overlap so providing a powerful in¬ 
strument of calibration. The calibration procedure is described in all details; it is self consistent 
because relies on two independent measurements of the same quantity, being the RPC signals 
independently readout by the strips and by the analog pickup electrode (BP). No additional in¬ 
formation is needed. The error on the measurement of the number of particles (n) crossing 
the BP is essentially determined by the intrinsic fluctuation of the RPC signal, which is about 
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Figure 23: Gain (mV/strip) versus incidence angle 6 of the shower, with respect to the zenith, for the four typology of 
BPs. 


20% for a single streamer at the operational conditions of the experiment, while decreases as 
1/ V(n)- The contribution to the error of the electronic calibration is less than 0.5%, while the 
electronics itself shows a quite good stability, or a seasonal variation of the response which is 
less than 1%. The mean signal of the single particle, the gain, is determined with an uncertainty 
of about (2.3%)stat + (3.5%)sys on the single run base, or a few hours; its mean time derivative is 
+(4.3%)/year which also represents the detector stability. More than 95% of the 3120 BP in the 
central carpet of ARGO-YBJ have been calibrated with the described procedure. The study of 
the trigger efficiency has confirmed the analog detector to be fully efficient above 100 TeV for 
all kind of primary cosmic rays. The physics performance of the detector has been evaluated by 
studying the number of particles hitting the BP at the shower core position (PMax). The spec¬ 
trum of this variable, as obtained by data analysis of showers with core in a fiducial area of the 
ARGO-YBJ carpet, has shown no systematic effects related to the use of different scales, that 
is a perfect matching in the overlapping regions of this distribution. An absolute comparison of 
the distribution with the Montecarlo expectation shows a very good agreement within the exper¬ 
imental uncertainties. The correlation of PMax with the energy, for different primaries, has been 
investigated, thus confirming that ARGO-YBJ can efficiently operate up to PeV energies beyond 
the knee of the primary cosmic ray spectrum. 

12. Conclusion 

As in the case with all air shower detectors, ARGO-YBJ samples the cascades induced deep 
in atmosphere by primary cosmic rays. Information about primary composition and energy spec¬ 
tra is inferred from measurements of the properties of the secondary air showers. The analog 
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Figure 24: Differential rate of PMax for events with core in a fiducial area of the carpet and 6 < 15°, showing a very 
good matching between different scales (see text). The results from a Montecarlo simulation are shown for comparison. 


readout of the ARGO-YBJ RPC detector, put in operation on December 2009, allows for a very 
detailed and precise measurement of the number of charged particles around the core of air 
showers induced by > 100 TeV primaries. This imaging is achieved by means of an array of 
3120 pixels ( called BigPad) of 1.7 m^ each, densely distributed over a 93% active area of about 
5770 m^ . The most important features of this device can be summarized as : 

• each pixel is instrumented to count the number n of charged particles with a resolution 
20%/ V(n) + 4.4% where the first term is related to the fluctuation of the streamer for¬ 
mation in the gas gap of the RPC , while the second term accounts for the uncertainty 
in the absolute calibration functions and control of the environmental parameters. As far 
the long-term stability of the detector, the m.i.p. signal shows a mean variation of of 
-1-4.3%/year; 

• the dynamic range of each pixel extends as far as particle densities of 2xl0^/m^ and more, 
allowing to measure the number of charged particles around the core of air shower from 
PeV primaries with an accuracy < 5%; 

• the performance of this device does not depends on the air shower incidence angle at least 
up to 50 degrees. 

In consequence, since the position of the detector at high altitude makes possible to sample 
ground particles near the shower maximum, thus reducing significantly the effects of fluctuations, 
one can reasonably predict accurate measurements around the knee of the primary cosmic ray 
spectrum with unprecedented resolution. Preliminary results presented at ICRC 2013 Il48]l as well 
as the study of the PMax distribution discussed in this paper confirm the excellent performance 
of the detector and of its analog readout. 
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Figure 25: The trigger efficiency of the analog system as function of the energy E of different primary cosmic rays 
(protons (blue squares). He (red triangles), Fe(black dots)). 
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